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The wave-front-shearing interferometer may be used to test any converging wave 
front regardless of whether or not it is symmetrical. A mathematical operation is described 
that permits complete analysis of the data. This operation yields values of the deviations 
of wave fronts under test from a close fitting sphere. The reference surface may be chosen 
statistically so that the results are the deviations from a best fitting sphere. 



Introduction 

The wave-front-shearing interferometer [1, 2, 3] * 
provides a comparison of a converging wave front 
with a sheared image of itself. Similarly, the 
wave-front-reversing interferometer [4] provides a 
comparison of either a converging or plane wave 
front with a reversed image of itself. A detailed 
description of a method of analyzing the data from 
these interferometers, except for the special case of 
revolution symmetry, has not been published. 
Consequently, the absolute shape of a completely 
unknown wave front could not be obtained. 

The method of analysis involves a mathematical 
operation. This operation is similar for the two 
aforementioned interferometers. Its description is 
more clearly portrayed in its application to the 
wave-front-shearing interferometer ( = WSI) than to 
the wave-front-reversing interferometer (=WRI), 
although the method evolved from a study of the 
latter. There are, however, sufficient differences in 
the two operations to warrant separate treatments. 
Accordingly, part 1, which follows, will describe 
the analysis of data from a WSI and part 2, to be 
published later, will describe the analysis of data 
from a WRI. Part 3 of this series will describe 
the absolute testing of optical flats. This method 
is not new but has not been described adequately 
in the literature. 

Parti 

1 . Analysis of WSI Fringes 

The WSI will yield unique solutions of a wave 
front. The assumption of symmetry (i.e., that it 
have revolution symmetry) is unnecessary. Any 
convergent wave front may be tested along any 
chosen diameter or along the line of its intersection 
with any plane that passes through the point of 
convergence. Several equally spaced reference 
points are chosen along this line. The line of ref- 
erence points must be parallel to the direction of 
shear and the shear must be an integral multiple 

i Figures in brackets indicate the literature references at the end of this paper. 



of the separation of the points. This arrangement 
and spacing of the reference points causes their 
images, in the region of overlapping wave fronts, to 
coincide in pairs as is shown in figure 1A. 

The separation of the reference points will be 
defined as unity for this discussion. If the magni- 
tude of the shear is unity also, one adjustment of 
the system (one set of fringes) is sufficient to eval- 
uate the deviations of the wave front at the chosen 
reference points from a statistically chosen sphere. 
For a shear that exceeds unity, more than one 
fringe pattern (i.e., two or more shear values) 
must be used to obtain a solution. 

In figure 1, P v (v=0, 1, 2 ... N) represents the 
(N-\-l) reference points in one of the images, W, 
of the sheared wave front and P' v the corresponding 
points in the other image, W. Let <5„ (see fig. IB) 
equal to the deviation of the wave front at P„, from 
the corresponding point, T„, on a reference circle, C, 
(to be chosen later), e equals the angle between the 
two images, C and C' ', of the reference circle at their 
point of intersection and /* equals the distance from 
the intersection point to P . The distance, ju? is 
positive if P is below (in fig. IB) the intersection 
of the two circles and negative if above. The devia- 
tion, <5„, is positive if P v is on the concave side of the 
circle, and negative if it lies on the convex side of 
the reference circle. The distance from P to P v (or 
from T to T v ), measured along the circle, is always 
positive and is represented by v. 

The separation, S v (distance from T v to jTJ_i), 
of the two images of the reference circle at any 
pair of points, P v and P f v -i, (assuming unit shear) 
may be obtained from figure 2. The centers of the 
two circles are located at d and C[. The distance 
OiC[ equals 2h, R equals the radii of the two circles, 
p equals the distance from E to T V) E bisects the line 
CiC[ and </> is the angle subtended by (p+/x) at E. 
Applying the law of cosines (from trigonometry) to 
triangles EC[T' V ^ and EC X T V we obtain 

R 2 =h 2 +(p+S v ) 2 -2h(p+S y ) cos (90-0) and 

R 2 =h 2 +p 2 -2hp cos (90+0). 
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Figure 1 A Two images of a wave front sheared laterally relative to itself. 

B Illustration of the two wave fronts relative to the images of a reference circle. 




Figure 2 Relationship between parameters and other variables. 

On taking differences and solving for S v we obtain, 

S„=2A sin <£. (1) 

This equation may be applied rigorously, but if the 
focal ratio (radius of curvature of the circle divided 
by the diameter of the wave front) exceeds six, no 



significant error is introduced by replacing sin </> by 
its approximate equivalent, (v-{-n)/R. Also, since the 
angle e will always be quite small, the value of 2 A is 
approximately equal to R-e and the value of S v in 
eq (1) becomes 

S v =(v+pl)€. (2) 

If we let Q iV represent the separation of the two 
wave fronts, W and W, at P V) it is seen from figure 
IB that 



v=l, 2, 3 



-S l , = 5 p _i J rQ v —e(v + n) 
. . N. 



} 



(3) 



All quantities that represent optical distances are 
given in units of the wavelength that is used. Ac- 
cordingly, these quantities should be multiplied by 
the wavelength to convert to standard units of 
length. The quantity Q v may be observed directly 
if white light is used to adjust the zero order to 
a known point, P r , and then, using monochromatic 
light of known wavelength, observe the order differ- 
ence between points P r and P v . It is inconvenient 
and unnecessary to use white light. Let Q r represent 
the unknown order of interference at an arbitrarily 
chosen point, P T , in the fringe pattern (fig. lA) and 
q v the difference in order (number of fringes) between 
this point and P v . The order of interference, Q V} at 
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P v is, therefore, equal to (Q r +Q>)- This introduces 
another unknown into the equations of observation; 
but this may be absorbed by replacing the product 
€jlc with a new parameter, r, such that 



en=(Q r —r). 



(4) 



With these values for Q v and ejj. replaced by their 
respective equals, eq (3) becomes 



8 v = 8 v _ 1 + q v +r— ve 
v=\,2 . . . N 



(5) 



Equation (5) represents N equations that contain 
(iV+3) unknowns: namely, (N+l) <5/s plus the 
two parameters e and r. We need three additional 
equations relating these (iV+3) unknowns if a solu- 
tion is to be possible. It appeared that one could 
use a second set of fringes formed with a shear of two 
units, as shown in figure 3A. This would yield 
(iV— 1) additional equations of observations with 
only two new unknowns, i.e., two more parameters, 
e 2 and r 2 , to be added to the above-mentioned 
(iV+3) unknowns. If N is 6 or larger, we would 
have as many equations as unknowns and a solution 
would seem to be possible, but it can be shown that 
these equations are not entirely independent. 

The <5/s do not have significance until the refer- 
ence circle, C, shown in figure IB, is defined. Any 
circle maybe defined by three conditions. In analyt- 
ical geometry, these might be the coordinates of the 
center and the radius of the circle. A circle might 
also be defined as the one that passes through three 
given points. The circle of reference may be fully 
defined by three equations of condition. When 
these are combined with the set of observation 
equations (eq (5)) data from a single set of fringes 
is sufficient for a complete solution. 

The values of the deviations may be obtained 
without knowing the value of R. However, when 
testing a parabolic mirror, for instance, the devia- 
tions represent the difference between the parabola 
and a sphere. If comparison is to be made between 
the measured deviations of the wave front and the 
computed deviations of a parabola from a circle, 
the radius of curvature of the circle must be known. 
The radius, R, can be measured independently of 
the interferometer. The ideal image point, E, 
(fig. 2) or the mean point of convergence of the wave 
front, can be located and its distance from any 
chosen point can be measured directly. The center 
of the reference circle will be located at the ideal 
image point if the reference circle is chosen as the 
circle that best fits the observed wave front. This 
may be done most precisely by the method of least 
squares. The value of R is the measured distance 
from the ideal image point (point of convergence) to 
any chosen reference point plus the computed devia- 
tion of this reference point from the reference circle. 

In the test of a concave mirror, R is approximately 
equal to the radius of the mirror. In the test of a 
lens (simple or compound system), R is the distance 
from the point of convergence (image point) to the 
back surface of the lens. 





Figure 3A Wave fronts sheared 2 units. 
B Wave fronts sheared 3 units. 

2. Specifications for the Reference Circle 

There are three relatively simple methods of 
defining the circle of reference. Each method re- 
quires three equations of condition. These con- 
ditions involve three well-known principles: namely, 
the method of coincidence, the method of averages, 
and the method of least squares. 

2.1. The Method of Coincidence 

We may define the reference circle, shown in 
figures IB and 2, by requiring it to pass through 
any three of the chosen reference points. Thus, we 
might require that it pass through the two end 
points, P Q and P N , plus some other point near the 
axis of the mirror (or lens). If the points, P and 
P N , are too close to a dubbed edge, a better fitting 
circle is obtained by requiring it to pass through the 
two points adjacent to P Q and P N , respectively. 
These two choices for reference circle are made by 
equating the corresponding 8 v 's to zero. Any three 
of the points can be chosen to define the circle and 
the results would be the deviations of the wave front 
at all chosen reference points from the chosen circle. 

For illustration we will choose the fringe pattern 
shown in figure lA. There are 8 reference points. 
We will arbitrarily require the reference circle to 
pass through points P , P 3 , and P 7 . This is done by 
requiring that 

8 =8,=8 7 = 0. (6) 

For this illustration, eq (5) represents 7 equations 
of observation. These 7 equations, when combined 
with the 3 condition equations, (6), form a set of 10 
linear equations with 10 unknowns; and all <5/s that 
are not evaluated in eq (6) can be computed in 
terms of known quantities (the #/s and v's). Table 1 
and figure 4 show the results obtained by this 
method. 

2.2. The Method of Averages 

The method of averages will usually yield a better 
fitting circle than the method of coincidence, because 
of the averaging of the 8 v 's. This method requires 
an assignment of weights. A reasonable correlation 
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Table 1. Deviations of a wave front from reference circles that 
are chosen as follows: Column 8, by the method of coincidence; 
column 4) by the method of averages; and column 5, by the 
method of least squares 



V 


Qv 


s, 




MC 


MA 


LS 







0.00 

.23 

.15 

.00 

-.06 

.06 
.21 
.00 


-0.20 

.13 

.11 

.00 

-.05 

.04 

.14 

-.16 


-0.10 


1 _. 

2 

3 

4 

5 

6 ._ 

7 


2.74 
2.10 
1.70 
1.45 

1.30 
1.00 
0.30 


.14 

.07 

-.07 

-.13 

-.01 

.15 

-.06 




Figure 4 Graphical representation of the deviations of a wave 
front from reference circles that are chosen by the method of 
coincidence (MC), the method of averages (MA), and the 
method of least squares (LS) . 

of image quality, formed by the measured wave front 
from an optical, image-forming system, with wave- 
front shape would normally require that the weights 
should be proportional to some function of the dis- 
tance from the reference points to the center of the 
circular aperture of the optical system. This would 
allow the assignment of equal weights to equal areas 
of the wave front. Equal weights will, however, be 
used in this discussion. 

In previous applications of the method of averages 
[5], the deviations of the observations from a curve, 
which is to be fitted to the observations by adjust- 
ment of parameters, are equated to functions of the 
parameters. These are then grouped into as many 
groups as there are parameters to be adjusted. The 
algebraic sum of the deviations in each group is 
equated to zero. When the functions of the pa- 
rameters (eq (5)) are substituted for the deviations, 
there are as many equations as parameters and the 
parameters can be evaluated. Any manner of 
grouping the deviations will yield a solution. A 
very poor and sometimes useless solution will be 
obtained if the points that are allocated to each 
group are scattered among those of the other two 
groups. This author prefers to avoid any inter- 
mingling of the points of one group with those of other 
groups. 



The parameters of the reference circle shown in 
figure 1A are €, r and any one of the deviations, <$„. 

The parameter e (see figs. IB and 2), together with 
R, determines the distance, h, from the center, C h 
of the reference circle to the point of convergence, E, 
of the wave front. The center of this circle is 
displaced from the ideal image point, E (mean point 
of convergence) for two reasons. One is that the 
axis of the mirror or lens cannot be adjusted to pass 
absolutely through the image of the source, and, 
secondly, the interference fringes would be too broad 
for reading when the wave front is nearly spherical. 
Consequently, the angle e is adjusted to be large 
enough to provide fringes of a most convenient 
width for accurate readings. 

The parameter, fi, together with R, determines the 
direction of the point C\ from E, relative to the axis 
of the mirror. The parameter r is related to and 
replaces /x through eq (4). Consequently, the 
parameters r and e define the center of the circle of 
reference relative to the wave front. 

Any one of the 5/s may be chosen as a parameter. 
By adding 8 V to the directly measurable distance 
from E to P v we have the radius R of the circle of 
reference. We will choose 8 for the third parameter 
because of simplification in the formulas that follow. 

This relationship of the three parameters to the 
coordinates of the center of the circle and its radius 
is somewhat vague, involved, and perhaps not 
clearly explained, but it is sufficient to understand 
that they correspond to the three parameters that 
define the circle of reference. 

If equal weights are to be assigned, there must be 
an equal number of points allocated to each of the 
three equations of condition. In the example chosen 
above (fig. 1A), the number of points are not inte- 
grally divisible by three. However, we may consider 
that each of the observation equations represents 
three identical observations. The total number of 
observation equations is then 3iV and division by 3 
is now possible. The three equations of condition 
are: 

350+3^ +2d 2 =8 2 +38 3 +38i+5 5 =28 5 +3dz+38 7 =0. (7) 

The set of equations represented by (5) may be 
replaced by the equivalent set 

d>=So+ib Zo+v • r-e S °-=jk (8) 

(7=1 <T = 

where j v is a function of the three parameters, <5 , 
€, and r. Any equation in set (8) is obtained by 
adding the first v equations in set (5) . The functions, 
jf„, may now be substituted for the corresponding 
5/s in (7). It will be noted that there is no f 
corresponding to 8 . This, however, does not present 
any difficulty because 8 is itself a parameter and 
may be left in eq (7) as is. 

The three equations of (7), after replacing <5„ by 
its equal, /„, contain the three parameters as un- 
knowns which may be evaluated. The computed 
values for the parameters are then substituted in 
eq (8) for evaluating the remaining <5/s. The results 
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for the chosen example, shown in figure lA, are 
given in table 1, column MA, and are also represented 
graphically in figure 4. 

The number of reference points, in the unit-shear 
method, is determined by the shear angle and di- 
ameter of the element being tested. For testing 
elements (mirrors or lenses) during the polishing 
and figuring operations, the optician is only inter- 
ested in departures from a desired figure. The 
easiest method of computing this is, therefore, the 
most practical. The shear may be adjusted to such 
a value that the resultant number of points is some 
multiple of three. The computation is a minimum 
when the number of points is the same in each group. 

2.3. The Method of Least Squares 

The method of evaluating a wave front by least 
squares is more elaborate than the others but would 
usually be used when highest precision is desired. 
In applying this method, the sum of the squared 



deviations 



(§0 



is to be minimized with respect 
required to define the 



to the three parameters 
reference sphere. 

This application of least squares differs from the 
usual method of least squares in that the observations 
are assumed to be free from error and the number 
of observations needs not exceed the number of 
unknowns whose values are sought (namely the 5's). 
In fact, a unique solution may be obtained even if 
the number of observations is one less than the 
number of chosen reference points. 

To minimize the sum of the squared deviations of 
the wave front from the reference sphere we require 
that 

N N 



X) 8l=$+y~*. /«=minimum. 



(9) 



This is effected by equating to zero the differentials 
of eq (9) with respect to r, e, and 8 0) respectively. On 
performing the differentiation and dividing through 
by 2, we obtain 



N bf N d/ N bf 



0. (10) 



bd £i J ' dr f^\ JV de 
The partial differentials, as obtained from eq (8), are 



d8 



1,(7=0,1, 2... N) 

a/ 



v =v, and > 
br de 



v-\ 

=-2 *, (f= 

<r=0 



= 1,2,3. ..iV) 



(ID 



On substituting these differentials into (10), we have 
for the example of 8 reference points (fig. lA), 



«o+«i+»a+«i+«4+««+*s+«7=0 

5x+2fi a +353+4fi4+55 6 +65«+75 7 =0 

52+3<5 3 +6<54+10<5 5 +15<5 6 +21<5 7 =0 



(12) 



On combining sets of eqs (8) and (12), we again 
have 10 equations with 10 unknowns. Using the 
same set of fringes (fig. lA) as was used above for 
illustration, we obtain for the <5„ 's the values given 
in column LS of table 1. These values are also 
plotted, along with those obtained by the other two 
methods, in figure 4. 

3. Sensitivity 

The sensitivity of the WSI varies with the absolute 
magnitude of the shear. If the shear is maintained 
at unity the sensitivity of the observations increases 
with a decrease in the number of points because the 
unit, as defined above, increases. The reduction of 
the number of points is, of course, limited. 

In testing surfaces that depart very greatly from 
a sphere, as for a large aperture and large aperture- 
to-focus-ratio parabolic mirror, it is desirable to 
reduce the sensitivity. Otherwise, such a mirror 
cannot be tested interferometrically with the WSI. 
As the number of points is increased, the absolute 
spacing of the points decreases. The smaller the 
si i car relative to aperture, the smaller can the optical 
path difference be made. By reducing the absolute 
magnitude of the shear the maximum optical path 
difference may be reduced at will and, consequently, 
the fringes may be adjusted to any desired width. 

The sensitivity may also be increased, without 
reducing the number of reference points, by using 
two or more shear values that are greater than 
unity (fig. 3). This may yield more equations than 
there are unknowns, but averaging of computed 
points is always acceptable. EacJi additional set of 
fringes, obtained with a different shear value, 
introduces two additional unknowns — new values 
for e and r. 

4. Symmetrical Surfaces 

Optical surfaces that are generated mechanically 
generally have revolution symmetry about the 
optical axis. If this is known and the axis is known 
to be at the center of the mirror (assuming circular 
apertures) then the deviations are known to be 
symmetrical because the center of the reference circle 
can be assumed to lie on the axis. Two other con- 
ditions are then sufficient to define the reference 
circle. The condition for revolution symmetry is 
defined by the specification 



Op fyv-y 



(13) 



and the condition equations, (6), (7) and (12), be- 
come, respectively, 

5 =8 3 =0, 



35o+35 1 + 25 2 = 65 3 +2$ 2 =0, 



and 



8 +8 1 +8 2 +5 3 =0 

21«$o+15<$ 1 +115 2 +95 3 



=»} 



(14) 

(15) 

(16) 
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The third equation in each of the sets (6), (7), and 
(11) becomes identical, respectively, to the first 
equation of the set. 

The assumption of symmetry (eq (13)) reduces 
the number of 5„'s by approximately ){N (exactly l / 2 N 
if N is an even number) and reduces the number of 
condition equations by one. If N is larger than 2, 
the number of equations exceeds the number of un- 
knowns, for a unit shear arrangement. Conse- 
quently, larger shears (2 or more units) may be used 
with increased sensitivity, as explained above. 
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